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Accretion disks in close binary systems are subject to a tidally driven 
parametric instability which leads to the growth of internal waves near the 
outer edges of such disks ( Goodman 1993| ). These waves are important in 
understanding the torque exerted on a disk by tidal forces and may play a role in 
the structure of the disk at small radii. Here we calculate the growth rate of this 
instability, including the effects of vertical structure and fluid compressibility. 
We find growth rates which are only slightly different from Goodman's original 
results, except that near the vertical resonance radius the growth rate can have 

O 
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vertical direction. Higher order modes, in the sense of increasing number of 
vertical nodes, have similar growth rates. Our results differ from a previous 
calculation along these lines by Lubow et al. (1993). The difference is mostly 
due to their neglect of radial gradients in the tidally distorted streamlines. 
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Introduction 



In a recent paper, Goodman (1993) pointed out that a tidally perturbed accretion 
disk would be subject to a parametric instability driven by the ellipticity of the distorted 
flow lines. This instability drives the creation of a standing wave pattern of low frequency 
waves with an amplitude that rises sharply near the outer edge of the disk. These waves are 
G-modes, or buoyancy waves, whose restoring force comes partly from vertical stratification 
of the fluid due to the vertical entropy gradient and partly from the radial stratification in 
the disk due to the steep radial gradient in specific angular momentum. Waves dominated 
by the former force are usually referred to as G-modes or (in a terrestrial context) internal 
waves. Waves dominated by the latter force are usually referred to as inertial waves. Here 
we will refer to them as internal waves. An overview of their properties in disks can be 
found in Vishniac and Diamond (1994, see also the references contained therein). 
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The generation of these waves is important for at least three reasons. First, the 
waves created in this manner carry negative angular momentum, relative to the angular 
momentum of the disk, and so their generation corresponds to a torque on the disk. This 
torque acts as a sink for the angular momentum transported outward in an accretion disk 
and will be important in determining its outer radius. Second, this process will lead to a 
strong inward flux of such waves in regions where the tidal forces are weak. The dissipation 
of these waves at small radii provides a weak, but potentially significant mechanism 
for angular momentum transport ( |Vishniac fc Diamond 1989) . Third, these waves can 



generate dynamo activity and angular momentum transport which, in an ionized disk, can 
dominate the angular momentum transport mediated directly by the waves ( |Vishniac, Jin, 



fc Diamond 199fl| , |Vishniac fc Diamond 1992] ). This last process will be insignificant if the 



Balbus-Hawley instability drives a dynamo with a growth rate comparable to fl QBalbus 



fc Hawley 1991| ) but this claim is inconsistent with phenomenological models of ionized 



accretion disks in binary systems (see, for example |Smak 1984a| , |Smak 1984b| , [Meyer fc" 



Meyer-Hoftmcister 1984 IHuang fc Wheeler 19891 , [Mincshige, fc Wood 19891 , |Mincshigc & 



Wheeler 19891 , |Cannizzo~1994l |Cannizzo, Chen, fc Livio 1995|) 



The physical basis of instability is closely related to a three dimensional instability 
of an incompressible fluid with elliptical flow lines ( |Pierrehumbert 1986| , [Bayly 1986| , 



Landman, fc Saffman 1987). An observer moving with the fluid in an accretion disk will see 



the tidal forces as inducing a set of time-dependent stresses whose comoving frequency, u>, 
is m(Q — Qf,), where m is the dimensionless azimuthal wavenumber of the tidal force, fl is 
the local orbital frequency in the disk, and fib is the binary orbital frequency. A realistic 
tidal force will have components at all values of m, but the strongest components will 
have small values of m. If we consider the m = 2 component then we see that an internal 
wave with ui = Q — Qi, and m = 1 can interact with this component of the tidally induced 
stress to produce a wave traveling in the opposite direction, with the same m and uj. In 
Goodman's paper he showed that this process could lead to exponential growth as each 
component of the standing wave amplifies the other. Since this process is inherently local 
no special boundary conditions are required. It is necessary that the waves be confined 
to the midplane of the disk, since the waves will need to maintain coherence over a radial 
distance comparable to their group velocity divided by the instability growth rate, but 
internal waves are naturally confined in this manner. 

The significance of this instability depends on its growth rate. A small growth rate 
not only implies weaker waves, but makes the instability vulnerable to suppression from 
dissipative effects, including turbulent dissipation from any other instabilities present. 
Goodman (1993) calculated the growth rate to leading order in the tidal forces for an 
adiabatic fluid in a disk without vertical structure assuming that the induced fluid motions 
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were incompressible. He found that at small r the growth rate was given by 



where M 2 is the mass of the secondary (mass-losing) star, and a is the semi-major axis 
of the binary orbit. This result is consistent with the growth rate found by Ryu and 
Goodman (1994), who used a different approach and who also carried out a two dimensional 
simulation to look at the nonlinear saturation of the instability. Goodman also found a 
singularity in the growth rate near the vertical resonance discovered by Lubow (1981). This 
resonance occurs at the radius where 



where 7 is the adiabatic index of the gas in the disk. This singularity is due to pumping of 
the waves by vertical motions rather than the ellipticity of the flow lines, but is a natural 
part of the parametric instability in accretion disks. 

More recently, Lubow, Pringle, and Kerswell (1993) repeated the growth rate 
calculation including the effects of an isothermal vertical structure. They found a reduced 
growth rate at all radii, by about an order of magnitude, which they ascribed to the 
stabilizing effects of the vertical stratification. However, their calculation differed in other 
respects as well. In particular, they neglected vertical motions, which eliminated the 
singularity, and they treated the flow lines as locally self-similar, which eliminated radial 
derivatives of the tidal distortions. They also suggested that the instability could be 
suppressed by a suitable choice of radial boundary conditions, but further work by Ryu et 
al. (1995) has shown that the instability is purely local and relatively insensitive to radial 
boundaries. 

In this paper we present an improved calculation of the linear growth rate of this 
instability, including the effects of vertical structure, compressibility , and variations in 
the local adiabatic index 7, while retaining tidally induced vertical motions and radial 
derivatives of flow lines. We show that a full calculation is reasonably consistent with 
Goodman's original estimate for the growth rate of the parametric instability in a realistic 
disk. However, there is one interesting difference. Vertical stratification greatly increases the 
amplitude and width of the peak in the growth rate associated with the vertical resonance. 

In §2 we derive the necessary formulae for calculating the linear growth rate. In §3 
we present our results for disks with different vertical structures and adiabatic indices and 
compare our results with previous work. In §4 we discuss the implications of this work and 
some directions for future research. 
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2. The Growth Rate Formalism 

The response of an accretion disk to tidal perturbations is an example of the more 
general situation where one has a set of modes of a system and one perturbs the evolution 
operator. In this case the basic modes are the m — ±1 standing internal waves of an 
accretion disk and the perturbation is the change in the evolution operator for these waves 
due to the tidal distortions. If a frequency eigenstate of the system is represented by some 
state vector A* then its unperturbed evolution is given by 

d t Ai = LijAj = iuAj, (3) 

where we have used implicit summation and uu is the comoving frequency, i.e. we treat the 
mode evolution from the point of view of an observer moving with the fluid. It is useful to 
define an adjoint state vector and corresponding eigenmodes, Aj, using the adjoint of L^, 
i.e. 

dtAi = -LljAj = icuAj. (4) 

When %L is a hermitian operator the two sets of eigenvectors are identical. We can define 
an inner product, (AA), in the usual way as 

(AA) = J A*A t dV, (5) 

where the integral is over the volume of the system. 

When the system is perturbed so that L — * L + AL then the eigenf unctions are 
perturbed as well. We can invoke the normalization constraint, i.e. 

(AA) = ((A + A)(A + A)), (6) 

to show that 

(A + A\L + AL\A + A) = i(Q + 5uj)(AA) = iuj(AA) + (A\AL\A), (7) 

or 

(A\AL\A) 

(AA) y ' 

In order to solve for the growth rate we need to find the unperturbed eigenvectors A and A 
and calculate the matrix elements due to tidal perturbations of the disk. 



We will assume that we are dealing with standing waves whose radial structure can 
be approximated by a plane wave of the form cos(/c r r + 0), where is an arbitrary phase 
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function, and whose azimuthal wavenumber is ±l/r. The latter is a requirement for an 
exact resonance with the m = 2 component of the tidal force. More generally we do not 
need an exact resonance to drive the instability, but we will restrict ourselves to that case 
in this paper. Waves with a low ker = m of order unity will show secular radial variations 
on scales ~ r/m flVishniac fc Diamond 1989|) . As long as we restrict ourselves to k r ^> kg, 
as we do here, then we can ignore this effect. 

We will restrict ourselves to adiabatic perturbations of a disk with arbitrary vertical 
structure. The evolution of radial plane wave perturbations of the disk, with radial 
wavenumber k r , are described by the following equation. 
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where a prime denotes a derivative in the z direction, S = P 1 ^ / p, c s (z) is the sound speed, 
p(z) is the density, P(z) is the pressure, 5 is the fractional density perturbation, and x is 
the pressure perturbation divided by P 1 ^. 

If we take the perturbation to have a well-defined comoving frequency, consistent with 
our approximation of the perturbation as a plane wave, we get the following second order 
equation in z for x 
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where A^ 2 is the square of the Brunt- Vaisala frequency, defined by 

N 2 = zQ{lnSy. 



X 



0. 



(10) 



'IV. 



This equation has an apparent singularity when uj 2 = N 2 , but closer examination shows 
that there are always two linearly independent and regular solutions passing through this 
point. The only boundary condition is the requirement that x vanish as \z\ — > oo. The 
other variables can be found in terms of x using equation @. Since terms of order mj (rk r 
are negligible, they are 
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There will be an infinite number of solutions that satisfy this boundary condition, just 
as there are an infinite number of harmonics of the resonant mode of a box. However, 
modes with many vertical nodes, or alternatively modes with large k z , will be more easily 
saturated due to nonlinear processes. Here we will consider only a few of the lowest order 
modes. In order to study the parametric instability we need to look at the standing wave 
pattern formed by the superposition of outgoing and ingoing waves. For a given uo and k r 
this involves an overall phase factor of the form 



X oc cos(fc r r + 0) cos{ujt + 9 + ip), 



(16) 



where <fi and r/> are phase constants and 9 is the azimuthal coordinate. The angle <fi is 
irrelevant and may vary from one annular region to another. (The instability only requires 
a limited amount of radial coherence.) Phase factors for the other components of the state 
vector can be derived from equations fli"2|) through (15). Given a model for the vertical 
structure of a disk the vertical eigenfunctions can be calculated from these equations in a 
straightforward manner. 

The evolution of the adjoint state vector is given by 
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where the tilde variables are defined by this equation. Note that we have used the relation 

(7-1) 
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in deriving equation fliTf ). The vertical structure equation, analogous to equation 
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Equation (|19D has an apparent singularity at G = 0. However, as before, one can show that 
there exists two linearly independent and well-behaved solutions that pass through this 
point. The other adjoint variables are given by 
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These equations can be solved, once again invoking x ~~ * as \z\ — > oo. Alternatively, 
we can notice that there is a close relationship between the state vector and the adjoint 
state vector. The fact that (AA) is a constant suggests that we should be able to equate the 
state norm with some conserved state property, e.g. the energy. In fact, if we compare v r 
and (1/2) pv r then we can show that equations fllQ|) , fll2|) , (|i~9|) , and ( ^l]) imply that the two 
variables satisfy the same second order differential equation in z. The boundary condition 
at \z\ — > oo does not distinguish between solutions that differ only by a complex factor, but 
the phase and amplitude of the adjoint state vector is undefined in any case. We conclude 
that we can take v r = (1/2) pv r . With this choice, we can use equations ( fl"2|) through ( |i~5| ) 
and (|1|) through (H) to show that 
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The meaning of this choice is clearer if we note that equation (|25D implies that 
(AA) = -p (\v r \ 2 + 4\ Vg \* + \ Vs f) + l^-JL. - -/r.0-' A; ! 
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Furthermore 




-±pzQ 2 Az* {-Azd z lnS) 

\ P N 2 \Az\ 2 . (28) 



Therefore, 



-p(KI 2 + 4N 2 + ki 2 ) + ^ ■ 2 . 



(AA) = t p ( \ Vr \ 2 + 4 |^| 2 + M 2 ) + J£-L + ^JV*|A*| a , (29) 



which is to say that using this convention the norm of the state is the locally measured 
energy density. Goodman (1993) used the local energy density in his calculation of the 
perturbation growth rate, but, consistent with his neglect of compressibility and vertical 
structure, left out the last two terms. 

It remains to find the operator AL and calculate the matrix elements (A|AL|A). 
We will start with the lagrangian displacements produced by the tidal forces. We 
follow Goodman in using first-order epicyclic theory to calculate the radial and angular 
displacements. The radial and azimuthal displacements are 

1 fd®2 2ft \ 

and 

-1 /2ft d$ 2 4cj 2 + 3ft 2 \ , n 
2(4ur — \l z ) \ru) dr r z uj z J 

where u = ft — ft&, 9 is measured relative to the axis running between the two stars, and $2 
is the azimuthal m = 2 Fourier component of the tidal potential. It is given by QBrouwer fc 
Ulemence 196 1| ) 



$ 2 = -^6^ 1/2 V/a), (32) 

(X 

where q is the mass ratio for the system, Mi is the primary stellar mass, and b^ 2 (r/a) is a 
Laplace coefficient. The vertical displacement follows from examining the response of the 
disk to long radial wavelength perturbations to the local vertical gravity. This gives 

zcos{2cu + 29) ( 1 \ 1 d (d<& 2 \ 4 3^ P _ 2 



^ _ zcos(2gJ-2^ / _ 2 /l g | 1 d ( d ® 2 ] 1 4 3, _ 

z Auj 2 — (7 + l)ft 2 \ \r r r r ) r dr \dr J r 2 r 

'(33) 

This differs from the expression in Goodman (1993) in two respects. First, the sign of the 
$2 terms are reversed. Second, there is an additional term proportional to £ r which comes 
from the change in vertical gravity induced by radial motions. Writing these displacements 
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as a set of coefficients of the angular functions cos(2co>/!; + 28) and sin(2u;t + 28) we find the 
tidally induced eulerian velocity perturbations are 



Av r = (Q£ e - 2u)£ r ) sin(2cjt + 28) 
Av e = f^Cr + 2&£o) cos(2cot + 28) 



(34) 
(35) 



and 



Av z = -2u£ z sin(2utf + 28). (36) 

In addition, we need to know the tidally induced perturbations to the pressure, density, and 
entropy. These are 
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P ' \ S p J 

Here we have used the assumption that the tidal perturbations are adiabatic. This 
assumption may not always be valid, but since we have already employed it in deriving the 
dispersion relation for the waves, we have not further compromised our results by using it 
here. In any case, this assumption should be reasonable as long as the cooling time for 
the disk is much greater than the orbital time. The radial and vertical derivatives of these 
expressions can be calculated in a straightforward manner. 

It is convenient to rewrite the components of the state vectors as coefficients of phase 
factors, following equation (fIB|). For both the state vector and its adjoint the phase factors 
have the form 

/ sin(/c r r + 0) sin(d)t + 8 + ip) \ 

— sin(fc r r + 0) cos(a;t + 8 + ip) 

- cos(k r r + 4>) sm(ujt + 8 + tp) . (40) 
cos(A; r r + 0) cos(a)t + 8 + ip) 

V cos(&vr + 0) cos(a>t + 8 + ip) J 
The definitions of the state vector and adjoint state vector components are then modified by 
dividing through by i for every factor of uj or k r . Note that when we are done constructing 



-10- 



our standing wave patterns there is no distinction left between A and A*. Finally, we note 
that matrix elements which involve products of the form cos(k r r) sin(fc r r) will not vanish 
exactly. Instead, 

(sm(k r r)F(r)cos(k r r)) = ~{k- x d r F{r)). (41) 

When the matrix element involves a factor of k r such terms will contribute to the total 
growth rate at the same order as other terms. This point is particularly significant since the 
radial gradients of the tidal perturbations can be much larger than the tidal perturbations 
divided by r. We note that in performing the necessary integration by parts we have 
neglected secular radial variations in the eigenmodes, which will arise from radial variations 
in density and pressure. 

We are now ready to write down expressions for the matrix elements which enter into 
equation (|8|). The tidally perturbed radial acceleration equation is 

d t v r = 2Qv e - )d r P - (Av r d r + ^d e + Av z d z + d r Av r ) v r + - ±d e Av r ) v e 

-v z d z Av r + tfSd rX +±d r AP. (42) 

Taking into account the radial phase factors given in equation (|16D we can use this to write 
the top row of AL as 

AL ri = (--d r Av r - Av z d z - —d e , - -d e Av r , 0, ^Sd r , ] . (43) 
V 2 r r r p J 

This gives a contribution to the growth rate, following equation (^ of 

(v r \AL\A) = ^f^v^ Av z d z v r - (2^4 - 2£&) Ve 

- x k r Ae. S y z 

= ^ IZo Vr U (ujd r i r - + (ffi - 2w) & - 2&f) + d z V r 2iu£ z 

_^(£ p (» + 2G>)+&(2fl-fi)) 
+ X k r S + d r i r + % + d z ^ z (1 + zd z \np)))dz. (44) 



Similarly, the tidally perturbed azimuthal acceleration equation is 
d t v e = -f v r - (d r Av g + ^) v r - (Av r d r + ^d 6 + Av z d z + ^ + 



t'9 



-v z d z Av e , (45) 
which implies a contribution to the growth rate of 

(vg\AL\A) = s *f& JZo v e U (d r Av e + + v e {-\d r Av r + ^ - ^) 
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+Av z d z vejdz 

= ^ J^oo v (v r (f <9 r £ r + 2ud r i e - + (2Cu - 3fi)&) 
+u« (-? drb + c^ r + (Jfi - 2Q) & - 2(Q + Q)Sfj 

-2uj£ z d z v^jdz. (46) 

The tidally perturbed vertical acceleration equation is 

8 t v z = -S(d z ln(l)+d z )x-zQ 2 S 

-v r d r Av z - ±<% Av z - v z d z Av z - (Av ■ v) v z 
+A P S + d z ) X + 5 -d z AP + 25^zQ 2 , (47) 

which implies a contribution to the growth rate of 

(v z \AL\A) = v z (v z (-d z Av z + \d r Av r - Am) - Av z d z v z 



= ^ roc V, (v, (Zdrto - C0d r Cr + &f - (ffl + 2w) & + 20^,) 

((7 - 2) ^ + 7 ^ - d£ z (2 + lnp + ^f)) S)dz. (48) 

The tidally perturbed pressure evolution equation is 

dtx = -§d r v r + ±(zn 2 v z - c 2 s d z v z ) 

+P~ 1 ^((-v-V)AP - 7APV • v- Av-V{P lh x) - 7V • {Av)P l ^ x y (49) 

which implies a contribution to the growth rate of 

(xl alia) = ^x{- c WKv r + 1 + ^O z ) v z 

((7 - I) - ^) + + d z ) + -yd,(Av,)) X )dz. 

= ^ He X(-| *r^«r - ^ " fce. (2 + ^ lnp + f )) „, 

+f¥4r + x(2fi>(^-7)fl.t+ 

(7 - i) (2fi 6 ^ - 6 (f + 4o)) + ft0 r & - 2w9 r e r )) - 2w^g)^ (50) 
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Finally, the tidally perturbed density evolution equation is 

d t 6 = -d r V r - V z d z hl(p) - d z V z - £f (d r V r + ^) - d r (ff) V r 

-^d e v e - v fd r (f) - fd z v z - ^d z Ap - 5(V • Av + Av z d z lnp). (51) 
which implies a contribution to the growth rate of 

(5\AL\A) = «^JZ o s(-4fk r v r + v M *ft + ^d z v z 

+ (±d r (rAv r ) - ^ + Av z d z lnp) 5 + Av^dz 
= ^ S-j(-^Kv r + v z (-d z { z (d z lnp - z{d z lnp) 2 + zdl lnp) + **d M lnp) 
+ T^f + 6 (H fi + 2u)^ + (Q-u)^ + §0 P £„ - w0 r £r - 2u;£A lnp) 

-2c^jW (52) 



In order to calculate the growth rate for a given vertical structure, we need to solve for 
the vertical eigenfunctions using equations fllOD through fllST) , and the adjoint eigenfunctions 



using equations ( |19| ) through ( P4|) or, more economically using equation fl25|), and substitute 
these results into equations (f44|) , (fiTf ), (H), (|50D, and (|52|) using equations (p7|), (|38|) and 
(j39p. Finally, we need to add up the matrix components and divide by the state norm given 
by 

(AA) = - [v r v r + vqVq + v z v z + XX + 86) dz. (53) 

4 J — oo 



3. Results 



We are now in a position to calculate the growth rates for the parametric tidal 
instability for a range of input parameters. Although our formalism allows us to take 
any vertical structure we will confine ourselves to presenting results only for adiabatic 
and isothermal vertical structures. These two cases should give us a sense of the range of 
possible results. We start with adiabatic disks. In this case the vertical structure is given by 
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and 



N 2 = 0. 



(58) 



The adiabatic index for the vertical structure is not necessarily the same index that appears 
in the perturbation equations. The latter describes the response of fluid elements on short 
time scales to changes in the ambient pressure, whereas the former is an attempt to model 
the vertical structure. In the preceding equations we have distinguished the two by using 
7 S when we mean the vertical structure parameter. This class of models has the peculiar 
property that the boundary occurs at \z\ = H rather than \z\ = oo. Consequently we have 
to modify our boundary conditions and take the solution which is least divergent near 
\z\ = H i.e. x g° es to a constant value near the vertical disk boundaries. An unfortunate 
consequence is that the fluid velocities are also constant, while the sound speed goes to 
zero. This is manifestly unphysical, but not too surprising. One expects the adiabatic 
approximation to break down far from the disk midplane. Fortunately the results of our 
calculation are weighted by the mode energy, and therefore by the local density, which 
drops to zero faster than the sound speed. It seems safe to assume that adding on an 
isothermal atmosphere will have little impact on our results as long as the column density 
in the atmosphere remains a small fraction of the total column density. 

In figure 1 we show the growth rate for the parametric instability in units of qVt^ for 
an adiabatic disk with 7 S = 7 = 5/3. as a function of r/a, where a is the semi-major axis 
of the binary orbit. In this figure we plot the results for the three lowest order modes of 
internal waves. The primary mode has one zero in x an d none in v z , i.e. the disk material 
bobs up and down as the wave passes by. In fact, in the limit where 7 S = 7 = 1, v z is not 
a function of z at all. (However, in this case v r does have a vertical mode. The material 
above and below the midplane are moving in opposite directions radially during each cycle.) 
We see that the second and third harmonics are almost indistinguishable in this plot, 
and the primary mode shows large differences only for radii close to, or greater than, the 
vertical resonance radius. The fundamental mode shows a stronger response to the vertical 
resonance and a significantly larger growth rate for larger radii. At small radial distances 
our results converge to Goodman's asymptotic formula for the higher order modes. The 
primary mode is weaker by about 23%. The vertical resonance peak is basically similar to 
the peak in his plot, but the zero in the growth rate occurs at a radius just outside the 
resonant radius, whereas in his calculation it was on the opposite side of the peak. The 
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source of this discrepancy is not clear. It is not due to our disagreement concerning the sign 
of the $ 2 terms in equation fl33|). 



In figure 2 we show the growth rate in units of qflb for the fundamental mode and for 
7 = 7 S = 5/3, 1.2 and 1.0. This last case is equivalent to taking an isothermal vertical 
structure and assuming isothermal perturbations. The two most significant effects of 
changing the adiabatic index are that the vertical resonance is pushed outwards and the 
peak around it is broadened. 

A more realistic, but still idealized, model is to assume an isothermal vertical structure 



for the disk. In this case we have P oc p, S oc p 



1-1/7 



and 



p = p exp I 2p I , (59) 

c 2 s = 7-, (60) 
P 



In figure 3 we show the growth rate in units of qQb for the first three internal wave 
modes for an isothermal disk with 7 = 5/3. We see that the response to the vertical 
resonance is much larger in this case, by a factor of about 5/2 for the fundamental mode. 
This carries over into a smaller response well beyond the radius of the vertical resonance. 
At much smaller radii the growth rates are similar to the adiabatic case, the only difference 
being that the primary mode is reduced from the higher order modes even more, by about 
36%. 

In figure 4 we show the growth rate in units of qflb for the fundamental mode in an 
isothermal disk with 7 = 5/3, 1.2 and 1.0. As before, lowering 7 broadens the vertical 
resonance peak and moves it outward. 

The modest differences between this calculation and Goodman (1993) for the adiabatic 
disk with 7 = 5/3 are explainable in terms of the additional effects included in this paper. 
On the other hand, our results are very different from Lubow et al. (1993). We have 
used a different calculational approach, and included the tidally induced vertical motions 
and pressure and density fluctuations. However, while these effects are quantitatively 
important, leaving them out does not produce a dramatic reduction in our growth rates, 
except near the vertical resonance. The most important difference seems to be our inclusion 
of radial derivatives in the tidal stresses. We can make a crude comparison between our 
two calculations by setting Ap = AP = £ 2 = 0, and taking d r ^ r = £ r /r and d r ^$ = £e/r. 



- 15 - 



In figure 5 we show the growth rates for the fundamental mode in an isothermal disk with 
7 = 5/3. The top curve shows the results of our complete calculation, also seen in figure 
3. The next curve shows the results of removing tidally induced vertical motions, and 
density and pressure fluctuations. The bottom curve shows the results adding to this the 
assumption that the streamlines are locally self-similar, i.e. the radial derivatives of £ are 
found by setting £ r and proportional to r in the interaction matrix. We note that at 
small r the bottom curve is more than a factor of 5 below the top one. This is a dramatic 
reduction, although this result still lies about 30% above the the reduction obtained by 
Lubow et al.. The remaining discrepancy is apparently due to the difficulty in entering 
precisely equivalent assumptions into our program, due to our somewhat different methods 
for evaluating the interaction matrix. We conclude that the bulk of the difference between 
the results of Goodman (1993) and Lubow et al. (1993) stems from the latter's neglect 
of radial derivatives of the streamline perturbations, not from their inclusion of vertical 
structure. 



4. Conclusions 



We have calculated the linear growth rate of a tidally induced parametric instability in 
accretion disks. Our results take into account only leading order terms in the tidal distortion 
of an accretion disk, but include compressibility, vertical structure, vertical motions, and 
radial derivatives of the tidal distortions. Our basic conclusion is that the original results 
of Goodman (1993), derived under a simpler set of assumptions, are basically correct. 
For a disk with adiabatic vertical structure we see only small deviations except near the 
vertical resonance, where we find the zero in the growth rate lies just outside the radius 
of the vertical resonance, rather than just inside it. This difference is unlikely to have any 
significant dynamical effects. For a disk with isothermal vertical structure we find a general 
enhancement of the growth rate near the vertical resonance, by a factor of about 2.5. This 
implies that the net torque on an accretion disk may depend significantly on its vertical 
structure. 

One lingering question regarding these results is the meaning of the infinite growth 
rate for the instability in the neighborhood of the vertical resonance. This singularity is a 
consequence of the divergence in £ 2 at that radius. In reality there are two effects which 
limit this divergence. First, strong vertical motions within a narrow annulus will lead to the 
emission of sound waves flLubow 1981|). First, equation (0) was derived neglecting radial 
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gradients. If we define an effective radial wavenumber for the growth rate by 

K r = d r \na(r), (62) 

then the tidal pumping of the vertical resonance will saturate due to the emission of sound 
waves from the resonant region when 

K r V group = « a{r). (63) 

Near the vertical resonance the growth rate is approximately 

r \ nl 



° (r)xA \Xr) q W (64) 

where Ar is the distance to the resonant radius and A is a function of 7 and the vertical 
disk structure. We see that the half-width of the annulus within which damping by acoustic 
emission is important is 

Ar ps ~ (65) 

2(fi - Q b ) qA{-f)n 2 b r qrA' 1 ' 

where H ~ c s /f2 is the disk thickness and we have ignored the distinction between Qj, and Q 
in the last (highly approximate) expression. Within this annulus the effective growth rate is 



2r 2 (n -n b ) ( nl 



r 



2 



viim^ "' V ^ 2 [Aq^j ~( q A) 2 Q b ^-) . (66) 

We see that as (H/r) — > acoustic emission becomes ineffective at removing the effects of 
this singularity. 

The second damping mechanism is due to nonlinear broadening of the vertical 
resonance. The resonance width will increase as the advective terms in the acceleration 
equation become important. In physical terms, the strong vertical oscillations of the disk 
will excite higher order harmonics, which will interfere with the resonant excitation of the 
fundamental mode. This is a purely one dimensional effect and does not depend on radial 
gradients in the tidal forcing. We can approximate its effects by inserting a term of the 
form (7o4u) 2 (<9 2 £ 2 ) 2 into the denominator of equation ( |3~3D where Co is a constant of order 
unity. This implies a limit on d z ^ z of order 




{d z t z ) Um ~ \ q-^ \ . (67) 
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Since the denominator of the leading factor in equation fl33|) goes as — 12Q b uo(Ar/r) 
this implies that this damping mechanism will be effective within an annular width of 
approximately 



2\ 2/3 



*»'m lBA *-~ r m[*3) • (68) 

The instability growth rate within this annulus will be limited to 

Comparing equations fl6"5|) and ( |68|) we note that acoustic emission will tend to be the 
most important damping process when the disk is not too thin, or when the tidal instability 
is driven only weakly There is no general rule as to which will be more important in 
accretion disks in compact binary systems, and it is possible that for some systems acoustic 
emission may dominate during outbursts while nonlinear broadening is important during 
quiescence. In either case the ratio of the torque exerted within the annulus dominated 
by the vertical resonance to the torque exerted on the rest of the disk at a similar radius 
can be estimated by comparing the cube of the growth rate weighted by the annulus width 
flGoodman 19931) . We have 



a(r) 3 r V 15 J \ArJ ' ^ 

For 7 = 5/3 A m 4.0 for an isothermal disk and 1.7 for an adiabatic disk. This suggests 
that the vertical resonance can dominate the torque on a disk and that accretion disks will 
have a tendency to have an outer radius that lies within this annulus. This tendency should 
be particularly noticeable when the disk structure is isothermal and/or 7 is relatively small 
(so that Q b /Q is relatively large). However, apart from the vertical resonance a(r) rises 
sharply with r and the torque per area, which goes as H 2 a(r) 3 rises even more sharply. 
This makes it possible for the dominant contribution to the torque on a disk to come from 
the outer edge of the disk provided that this lies a significantly greater radius than the 
vertical resonance. 

In this paper we have assumed that the orbital motion of gas within the disk deviates 
from perfect circles only to the minimal extent required by tidal forces. For most systems 
this is a reasonable assumption, but there are particularly cases where it is likely to 
fail. For example, SU Uma stars exhibit superoutbursts which are probably due to the 
expansion of the outer edge of the disk through the 3 : 1 orbital resonance ( Whitchurst 



|1988|) . Lubow (1991a, 1991b) has shown that at this radius the disk is vulnerable to an 
elliptical instability. Under these conditions the tidal parametric instability discussed here 
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can exist as a secondary instability, implying a dramatic rise in the generation of internal 
waves during a superoutburst. 
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Fig. 1. — Growth rates for the parametric instability in units of qQb as a function of r/a 
for an adiabatic accretion disk with 7 = 5/3. The solid line is for the minimal vertical 
wavenumber. The dotted and dashed lines are for the second and third wavenumbers. 




Fig. 2. — Growth rates for the minimal vertical wavenumber mode in an adiabatic disk for 
7 = 5/3 (solid line), 7 = 1.2 (dotted line), and 7 = 1.0 (dashed line). 




Fig. 3. — Growth rates for an accretion disk with an isothermal vertical structure and 
7 = 5/3. The solid line is for the minimal vertical wavenumber. The dotted and dashed 
lines are for the second and third wavenumbers. 




Fig. 4. — Growth rates for the minimal vertical wavenumber mode in an isothermal disk for 
7 = 5/3 (solid line), 7 = 1.2 (dotted line), and 7 = 1.0 (dashed line). 
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Fig. 5. — Growth rates for the minimal vertical wavenumber mode in an isothermal disk 
with 7 = 5/3. The solid line is full result. The dotted line shows the effect of suppressing 
vertical motion and tidal perturbations to density and pressure. The dashed line results 
from also assuming £ oc r. 



